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Understanding hadronic interactions is crucial for investigating the properties of unstable hadrons, 
since measuring physical quantities for unstable hadrons including the resonance mass and decay 
width requires simultaneous calculations of final scattering states. Recent studies of hadronic 
scatterings and decays are reviewed from this point of view. 

The nuceon-nucleon and multi-nucleon interactions are very important to understand the forma¬ 
tion of nucleus from the first principle of QCD. These interactions have been studied mainly by 
two methods, due originally to Liischer and to HALQCD. The results obtained from the two 
methods are compared in three channels, 1 = 2 two-pion, H-dibaryon, and two-nucleon channels. 
So far the results from the two methods for the two-nucleon channels are different even at the 
level of the presence or absence of bound states. We then discuss possible uncertainties in each 
method. 
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1. Introduction 

The ultimate goal of lattice QCD is to understand the properties of all hadrons quantitatively 
by first principle calculations. For stable hadrons, recent lattice QCD calculations have reached the 
level of a precise measurement of physical quantities such as mass and decay constant, thus touch¬ 
ing the goal. On the other hand, we are far from this goal for unstable hadrons. The calculations 
and analyses involved are much more difficult and complex, because one has to simultaneously 
treat final scattering states to measure physical quantities of unstable hadrons. To calculate phys¬ 
ical quantities in the scattering and decay processes, it is important to understand the interactions 
between hadrons based on QCD. 

The most famous of hadronic interactions is the nuclear force, which binds nucleons into 
nuclei. Thus, another goal of the lattice QCD is to prove the formation of nucleus from the degrees 
of freedom of quarks and gluons. 

For those goals, so far various studies have been carried out. In this report, recent studies on 
hadronic scattering, decay, and light nuclei are reviewed. Furthermore, it is reviewed consistency 
and inconsistency of results obtained from Luscher’s method and HALQCD methods in / = 2 two- 
pion, H-dibaryon, and two-nucleon channels. In the two-nucleon channels, possible systematic 
uncertainties in each method are discussed. 

The organization of this report is as follows. In the next section, Luscher’s method is briefly 
reviewed, which allows one to evaluate the scattering phase shift of two particles from finite volume 
calculations. In Secs. 3 and 4 some recent studies for the scattering length and physical quantities 
in decay processes are presented, respectively. Sec. 5 presents a compaiison of results from the 
Luscher’s method and HALQCD method, the latter having been recently proposed to calculate the 
potential between two hadrons, and a discussion on the uncertainties in each method. In Sec. 6 we 
review recent calculations of light nuclei and related studies. Conclusions are given in Sec 7. 


2. Ltischer’s finite volume formula 


Luscher’s formula [1, 2] in the original form enables to evaluate the scattering phase shift 8{p) 
of a two-particle system in the center-of-mass frame from the energy of two-particle state in finite 
volume, e.g., in the S-wave case. 


tan5(p) 


Zoo(l;^^) ’ 


(2.1) 


where q = Lp/{2n) and 


( 2 . 2 ) 


The relative momentum between the two particles, p, is determined from the two-particle energy, 
= 4(m^ -|- p^) with m being mass of the particle. An important assumption is that the interaction 
of the two particles does not depend on volume, and that the interaction is well localized, i.e., there 
is a region such that F(r) « 0 for r < L/2 with L being the linear spatial extent. 

An extension of Luscher’s formula to moving frame, where the total momentum of two par¬ 
ticles is non-zero, was first proposed by Rummukainen and Gottlieb [3]. Their result was con¬ 
firmed by different derivations [4, 5], and extensions in moving frame with different total mo- 


2 






Hadronic Interactions 


Takeshi Yamazaki 


menta, and several irreducible representations have also been derived [6, 7]. In the recent pa¬ 
pers [8, 9, 10, 11, 12] extensions of the formula in moving frames to the case of particles with 
different masses are derived, in which case mixings occur between even and odd partial waves. 
Those extensions are also reviewed by other plenary speakers [13, 14]. 

Luscher’s formula encompasses the case of a bound state as a pure imaginary solution for the 
relative momenta p, and indicates that the size dependence of the energy shift of the ground two- 
particle state relative to the free two-particle state is exponentially small [15, 16]. With this result 
as a background, study of size dependence of the energy shift has been employed in a number of 
recent papers to examine bound state formation in multi-hadron systems [17, 18, 19, 20, 21]. We 
shall call the methodology in these studies generally as Luscher’s method. 

3. Scattering length 

The scattering length qq is the simplest physical quantity in an S-wave scattering process, 
which is defined from the phase shift 5 (p) by qq = linip^o [tan ^ (p)/p] in the zero relative momen¬ 
tum limit. 

3.1 1 = 2 two-pion channel 

The 1 = 2 two-pion scattering is a very simple process involving no bound states or resonances 
in the low energy region. It is also the simplest from lattice calculation point of view, where only 
the direct and cross diagrams are necessary for the two-pion correlation function. While there have 
been various quenched calculations of ao in this channel, only dynamical calculations are reviewed 
in this report. 

The left panel in Fig. 1 shows the results for aoinji from dynamical calculations [22, 23, 24, 
25, 26, 27, 28, 29, 30] together with two experimental results by E865 [31] and NA48/2 [32] 
plotted with burst symbols. The calculations by Hadron Spectrum [33, 7] and NPLQCD [34] 
Collaborations mNf = 2 + \ full QCD are not plotted here, because the pion decay constant has 
not been calculated in the simulation parameter. The lattice results are compared with LO (leading 
order) ChPT (chiral perturbation theory) formula, a^m^i = — ni|/(87r/|), drawn by the dashed 
curve, where is the pion decay constant with the normalization of = 132 MeV at the physical 
pion mass. Most of lattice results agree with the formula. In the small ntj^jfji region, however, 
partially quenched staggered (open square) [28] and nonperturbatively improved Wilson (closed 
circle) [30] results significantly deviate from LO ChPT. 

The reason for the deviation of the staggered calculation [28] is not clear. The lightest data in 
Ref. [28] (open square) should agree with the one at the same m-ji !in Ref. [29] (closed square), 
because their actions and simulation parameters are the same. On the other hand, the deviation of 
the improved Wilson calculation is probably due to chiral symmetry breaking effect of the Wilson 
quark action at non-zero lattice spacings [30]. 

The lattice results for oq at the physical m^i are shown in the right panel of Fig. 1 . For the 
chiral extrapolation of qq, NLO (next-to-leading order) ChPT formula is employed in almost all 
calculations [23, 26, 28, 29]. On the other hand. Refs. [24, 30] use variations of NLO ChPT 
formula including systematic errors coming from finite lattice spacings. The value of Ref. [34] is 
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estimated using the low energy constants extracted from 8{p) atniTi = 390 MeV with NLO ChPT 
formula for the scattering amplitude, while they use the value of fji calculated with different action. 

The lattice results are compared with the phenomenological determination, marked “CGL”, 
obtained with NNLO ChPT [35] and the experimental results [32, 31] with and without ChPT 
constraint. Since in ChPT the scattering lengths for the 7 = 0 and 2 channels share a common low 
energy constant, the error becomes smaller in the case with the constraint. Some lattice results have 
smaller error, even if including the systematic error, than the phenomenological determination and 
the experimental results. Therefore, qq calculation in the 7 = 2 two-pion channel enters precision 
measurement era, where lattice result will be confirmed by experiment. In this conference, prelim¬ 
inary result of 77/ = 2 -|- 1 -|- 1 twisted mass quark calculation [36] is reported, which will calculate 
more precise value of qq at the physical lUji. 
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Figure 1: Left panel: aq in the 7 = 2 two-pion channel calculated in full lattice QCD as a function of fn, 
together with the prediction of LO ChPT (dashed curve). Right panel: values at physical m-^. Circle, up 
triangle, left triangle, down triangle, and square symbols denote results obtained with clover, domain wall, 
twisted mass, overlap, and improved staggered quark actions for valence quark, respectively. MA and PQ 
express mixed action with asqtad sea quark and partially quenched calculation. The experimental result by 
E865 [31] and NA48/2 [32], and the phenomenological determination using NNLO ChPT [35] (denoted by 
CGL in the right panel) are also shown. 


3.2 7 = l/lKn channel 

The S-wave scattering in the 7=1/2 Kn channel is more complicated than the 7 = 2 two- 
pion channel; it is in fact much more interesting because of the expected presence of fc meson 
from phenomenology in this channel. In the region of heavy quark mass, fC meson is indeed a 
bound state and so has to be taken into account in the ao calculation of Kn scattering [30]. In 
the case, one needs to calculate the first excited state in the system to evaluate qq- Furthermore, 
for the calculation of Kn correlation function, we need to include the rectangular diagram, for 
which specific calculalion melhods are required, for example Ihe slochaslic LapH melhod [37] or 
Ihe Irapezoid diagram calculalion melhod [38]. 

Due lo Ihese difficulfies, Ihere have been only four sludies in Ihe 7=1/2 Kn channel; 77/ = 0 
calculalion by Nagala et al. [39], partially quenched 77/ = 2 -|- 1 by Fu [40], 77/ = 2 by Lang et 
al. [41], and 77/ = 2 -|- 1 by PACS-CS Collaboralion [30]. The left panel of Fig. 2 shows Ihe resull 
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of the four studies for aQm„, together with phenomenological determination [42]. The lattice data 
do not show a common dependence in contrast to the situation with the 1 = 2 two-pion channel 
in Fig. 1. This is presumably due to the systematic error inherent in each calculations, e.g., a small 
spatial size of ~ 2 fm in the calculation of Ref. [41] and chiral symmetry breaking effect at the 
lightest data point of Ref. [30]. 

Similar to the calculations with the 1 = 2 two-pion channel, aoniji at the physical quark mass is 
estimated using NLO ChPT formula in Ref. [40] and its extension with chiral symmetry breaking 
term in Ref. [30]. These results are plotted in the right panel of Fig. 2. The panel also presents an 
indirect estimation by NPLQCD Collaboration [43], which uses the result of ao in the I = 2/2 Kji 
channel, since qq for / = 1/2 and 7 = 3/2 Kn channels share the common low energy constants in 
NLO ChPT formulae. 

Lattice results are in rough accord with the phenomenological determination. While this is en¬ 
couraging, in order to make more precise comparison with experiment, further direct calculations 
would be desired. In this conference, as shown in Fig. 2, RBC-UKQCD Collaboration [44] reports 
preliminary result of direct calculation of aq at the physical quark masses, where chiral extrapola¬ 
tion is not needed, in a large spatial extent of 5.5 fm. Furthermore, after the conference. Hadron 
Spectrum Collaboration presents the result at = 390 MeV [45]. 
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Figure 2: Left panel: aontj^ as a function of in the 7=1/2 Kn channel. Right panel: aoirijc at physical 
m„:. Phenomenological determination by Biittiker et al. [42], and an indirect determination by NPLQCD 
Collaboration [43] are also shown. Burst symbol in left panel is shifted for clarify. 


4. Resonances 

It is possible to obtain the mass and decay width of resonances from the phase shift 8{p) of 
its final scattering state, at least in two-particle decay processes. For the p meson decay into two 
pions, for example, the p meson mass nip and an effective pnn coupling gj^p can be extracted 
through the Breit-Wigner form. 


^cot5(p) 


^{mp-s), 

&nKp 


(4.1) 
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where ^/s is the center-of-mass energy of the two pions. The decay width Fp is then obtained by 


rp = 


Pp SpTCTC 
iffip 6n 


(4.2) 


where pp is the relative momentum of the two pions at the resonance point {s = m^p), =m^/4 — 

m\. Therefore, when 8{p) is calculated by Luscher’s method or its extensions for at least two 
different values of p, nip and gj^ip (or Fp) can be determined using eq.(4.1). 

In order to obtain the phase shift for several different values of p, we might need to extract 
the excited state energy corresponding to higher scattering states. For this purpose, the variational 
method proposed by Liischer and Wolff [46] is often used in lattice calculations. 


4.1 7=1 two-pion channel 

The decay process of the p meson is an ideal resonance process for an initial study of scattering 
and decay, since the spectrum consists only of the p meson and its final two-pion scattering state 
from the threshold (Im^) to the vicinity of the resonance energy region. Computationally, however, 
this decay process is much harder than the ao calculation of the 7 = 2 two-pion channel. One of 
the difficulties is the necessity of the rectangular diagram, as in the 7=1/2 Kn channel, with a 
non-zero relative momentum of the two pions. Another difficulty is that, since the decay proceeds 
through the P-wave, the pion has to be sufficiently light and the spatial size L sufficiently large so 
as to satisfy the decay condition, 

ml>A{inl + {lTi/Lf). (4.3) 

This condition can be relaxed if one uses the moving frame [3]. 

The first calculation of the p meson decay was carried out in A/ = 2 QCD by CP-PACS 
Collaboration [47] using the total momentum |P| = 1 frame. After this work, five papers [48, 
49, 38, 50, 51] have been published. Figure 3 show 8{p) calculafed wifh nijt = 390 MeV by 
Hadron Specfrum Collaboration [51]. There are 29 dafa poinfs obfained af differenf p’s using 
several moving frames and irreducible represenfafions on fhree differenf volumes. From fhe Breif- 
Wigner fif as in eq. (4.1), fhey exfracf fhe resonance mass mn = lUp and fhe coupling g^j^p as shown 
in fhe figure. 

Figure 4 presenfs a compilation of laffice resulfs for nip (leff panel) and gj^^p (righf panel) in 
Refs. [47, 48, 49, 38, 50, 51] exfracfed from fhe Breif-Wigner form fif. The experimenfal values 
are also ploffed. For nip resulf, fhe Nf = 2 resulfs by ETM Collaborafion are larger fhan fhe ofhers. 
This discrepancy could be explained by a sysfemafic error from scale seffing. If we choose Sommer 
scale ro for fhe scale seffing in all fhe resulfs, discrepancy becomes smaller, for example, fhe resulfs 
for ETM and PACS-CS Collaborations differ by less fhan 7% [38]. 

The laffice resulfs for gjtjcp, presenfed in fhe righf panel of Pig. 4, show fhaf if has only weak 
dependence on fhe pion mass, and fhaf ifs value is consisfenf wifh experimenf albeif errors are sfill 
generally large. We should nofe fhaf fhe resulf by Hadron Specfrum Collaborafion [5 1 ] wifh a small 
error does nol quife agree wifh experimenf, and so is fhaf by Pang et al. [49] albeif wifh a beffer 
agreemenf. The discrepancies from fhe experimenf are nof decreased by choosing vq fo sef fhe 
scale. In fhis conference, fwo preliminary resulfs by Pahy et al. [52] and BMW Collaborafion [53] 
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are reported as shown in Fig. 4. Fahy et ah [52] obtain the result at = 240 MeV, and BMW 
Collaboration obtain the results in = 134-300 MeV [53]. 



Figure 3; / = 1 tttt phase shift calculated by Hadron Spectrum Collaboration [51]. 
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Figure 4: Resonance mass trip (left panel) and the coupling gj^j^p (right panel) obtained in various calcula¬ 
tions. Open circle is result with lattice dispersion relation [47]. 


4.2 Other decay channels 

In the P-wave 7=1/2 channel, the K* meson decays to the Kk final state. This is an ideal 
decay process as in the p meson decay; there are only Kk scattering and K* state in the energy 
region below the resonance. There have been three lattice calculations of 5(p) in this channel, by 
Fu and Fu [54] in Nf = 2 + \ QCD, Prelovsek et al. [55] in 77/ = 2 QCD, and Hadron Spectrum 
Collaboration [56, 45] in 77/ = 2 + 1 QCD. The last one is reported in this conference, and reviewed 
by other plenary speakers [13, 14]. A difficulty of the lattice calculation in this channel is mixing 
of even and odd partial wave contributions in moving frame [8, 9, 10, 11, 12] due to mk- 
Some methods to solve the difficulty were proposed in Refs. [57, 55, 45]. 

In Ref. [54] it is assumed that S-wave contribution is negligible, so that the result contains sys¬ 
tematic errors from the mixing. Prelovsek et al. [55] obtain a value consistent with the experiment, 
gTiKK* = 5.7(1.6), at Mjc = 270 MeV. Hadron Spectrum Collaboration [56, 45] calculate 5{p) at 
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niji = 390 MeV, while the K* meson is a bound state at this value of mji. Hadron Spectrum Col¬ 
laboration also carry out coupled channel analyses in S- and D-wave channels to extract 8{p) for 
each of tiK and r]K scatterings and inelasticity. 

We should add that studies related to the decay channels, N* —)• Nk [58], a —)• nn [59, 60], 
and ao ^ TiTj [61] are reported at this conference. 

5. Comparisons of LUscher’s and HALQCD methods 

In this section, recent results obtained by Liischer’s method and HALQCD method are com¬ 
pared in three channels, 1 = 2 two-pion channel, H-dibaryon (AA) channel, and two-nucleon NN 
channel. After the comparison, possible uncertainties in each of the two methods are discussed. 

5.1 HALQCD method 

The HALQCD method [62, 63] extracts the potential between two particles, V (r), from Nambu- 
Bethe-Salpeter (NBS) wave function. The idea derives from the derivation of Liischer’s formula in 
quantum mechanics [2] and in quantum field theory [64], and also from calculations of the NBS 
wave function in two-dimensional scalar field theory [65] and 1 = 2 two-pion channel [66]. While 
these studies concentrated on the behavior of NBS wave function at large r, where V (r) « 0, the 
HALQCD method asserts that dynamics could be extracted from the smaller r region. 

The NBS wave function of the u-th NN state, ^n{r), in NN channel is calculated from four- 
point function of N operator, 

= L(0|A(Y+r,OA(Y,OAW(0)|0) + (5-1) 

X n 


where AA is a AA source operator, W„ is the u-th energy, Wn = 2yjm^ -|- p^, An = (AA, W„|AA|0), 
and 

Ur) = £(0|A(Y+r)A(x)|AA,W„). (5.2) 

X 

The spin index of the A operator is omitted for simplicity. The dots represent contributions from 
inelastic scattering states. A non-local potential Uir,f) is defined by acting with the Laplacian on 

Ur), 



Assuming that the higher order terms of the velocity expansion of U {r,f) are negligible, 

U{r,~r) ^V{f)5{f—~r), (5.4) 

a local potential V(j) is obtained from eq.(5.3). It should be noted that in quantum field theory 
the right-hand side of eq.(5.3) corresponds to Fourier transformation of off-shell scattering ampli¬ 
tude [66], which is related to 5{p) only on shell. 

One way to calculate V (r) is to use the ground state NBS wave function ^o(r) extracted from 
CNN{r,t) for large t. In the AA case, however, statistical noises increase exponentially at large t, 
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rendering an accurate calculation of CNN{r,t) difficult. It was proposed in Ref. [67] to employ the 
time dependence of Cj^i^ird) and rewrite eq. (5.3) as, 

/ 1 5 \ _ 

- ]CNNir,t) = V(r)Cm{r,t), (5.5) 

\mN ^niN dt^ dt J 

where C^Nird) = CA?Ar(r,f)exp(2mArt). A key assumption for this equation is that all (/>„(r)’s in 
CNN{r,t) yield a common V{f), in other words, eq.(5.4) is valid. With eq.(5.5), V{r) can be ob¬ 
tained from CNN{f^d) smaller t region as long as contributions of inelastic states are negligible. 

In the HALQCD method, thus, the potential V (r) obtained is fitted by a continuous function 
of r, like the Yukawa potential plus some hard core, which is then employed to calculate the phase 
shift S(p) by solving the Schrddinger equation. An assumption in calculating d(p) is that finite 
volume effects are negligible in V(r). It is also assumed that, while the potential V(r) may depend 
on the choice of the sink operator employed for C^Ni/d) iri eq-(5-l), such a dependence does not 
affect the phase shift 8{p) evaluated from the potential. 

In the left panel of Fig. 5 we present a recent result of HALQCD for the nuclear potential of the 
spin-singlet 'Sq channel calculated in A/^ = 2 -|- 1 QCD at ruj^ = 410-700 MeV [68]. The potential 
has the features similar to the phenomenological nuclear potentials with a repulsive core in the 
smaller r region and a one-pion exchange potential at large r. The depth of the lattice potential 
seems shallower than the phenomenological ones, which is considered to be caused by a heavy nin 
used in the calculation. The right panel of Fig. 5 presents the results of 8{p) at each calculated 
from V(r). 

In this conference, using the HALQCD method, an update of the Q.Q. potential [69] and a 
study of the quark mass dependence of the three-nucleon potential [70] are presented by HALQCD 
Collaboration. 




Figure 5: Recent result of nuclear potential (left panel) and 5 {p) obtained from the potential (right panel) 
in A/ = 2 -I- 1 QCD calculated by HALQCD Collaboration [68]. 


5.2 1 = 1 two-pion scattering 

Kurth et al. [71] calculate 8 (p) for the 1 = 2 two-pion channel by the two methods for Nf = 0 
QCD at = 940 MeV on four volumes from L = 1.8 to 5.5 fm. The stability of V (r) calculated by 
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the t-dependent method of eq.(5.5) is illustrated in the left panel of Fig. 6, where different colors 
correspond to different t ranges, t = 15^8(/?i), 24-48(/?2), and 33^8(/?3). The right panel of 
Fig. 6 compares 5{p) calculated from V{r) (solid lines) with that extracted by Liischer’s method 
(filled symbols). For this repulsive channel with a relatively simple shape of V(r), we observe an 
agreement of the two methods for small momenta. 



Figure 6: Potential in / = 2 tttt channel (left panel) and d (p) obtained from the potential (right panel) in 
quenched QCD calculated by Kurth et al. [71]. 


5.3 //-dibaryon channel 

HALQCD Collaboration [72, 73] examines the formation of H-dibaryon as a bound state in 
the AA channel mNf = 3 QCD. The result of V{r) at ruj^ in the range 470-1200 MeV obtained 
by the HALQCD method is shown in the left panel of Fig. 7, which exhibits an attractive core in 
contrast to the NN case. Fitting V(r) and solving the Schrddinger equation, a bound state is found 
for all values of m^i examined. 

The right panel of Fig. 7 shows the result from Liischer’s method for A/ = 3 at = 800 MeV 
by NPLQCD Collaboration [20]. The effective energy shift of the ground AA state from 2rn\ is 
clearly non-zero in the plateau region. After investigating the volume dependence of the energy 
shift, it is concluded that the ground state is a bound state. 

The binding energies calculated in the two methods are plotted in the left panel of Fig. 8 
where we also added the result of NPLQCD Collaboration for A/ = 2 + 1 [17]. Qualitatively, the 
two methods are in agreement in that both predict bound states for the range of m^i examined. In 
quantitative detail, they differ significantly at m^i ~ 800 MeV. 

In this conference. Green et al. report a preliminary result [74] in A/ = 2 QCD at mji = 1 GeV 
and 450 MeV. They employ the variational method [46] with local products of six quark operator for 
both source and sink time slices, and conclude that the ground state has higher energy than 2mA, as 
presented in the right panel of Fig. 8 for the case of m^i = 450 MeV, in contrast to NPLQCD result. 
A difference between NPLQCD’s and their calculations is sink operator. NPLQCD Collaboration 
uses two-baryon sink operator. After the conference, a preliminary result of calculation with a 
two-baryon operator is presented in other conference [75], and they observe smaller ground state 
energy than 2niA. 
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Figure 7: Potential calculated by HALQCD Collaboration [73] (left panel) and effective energy shift calcu¬ 
lated by NPLQCD Collaboration [20] (right panel) in A/ = 3 H-dibaryon channel. 
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Figure 8: Binding energy obtained from Liischer’s and HALQCD methods (left panel) and preliminary 
result of effective energy shift calculated by Green et al. [74] (right panel) in H-dibaryon channel. Woe and 
Vmax in the left panel express results in the infinite volume and on the largest volume, respectively. 
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5.4 Two-nucleon channels 

Two nucleons can be in either in the spin-triplet ^Si channel or in the spin-singlet ^Sq; in 
nature the former has a deuteron, while there are no bound states in the latter. 

Figure 9 summarizes all results obtained so far by Luscher’s method [76, 77, 18, 78, 20, 19, 21] 
for the energy shift of the ground NN state from Im^ in the spin-triplet ^Si (left panel) and spin- 
singlet *So (right panel) channels. Early studies [76,77] did not examine the volume dependence of 
the energy shift, so that it is not clear if the ground state is a scattering state or a bound state. Recent 
studies consciously examine this issue by working with multiple volumes [18, 78, 20, 19, 21]. 
Except for one result [78] (up triangle) which unfortunately suffers from large error, all results 
obtained by study of volume dependence have found that the ground states in both the channels are 
bound states. The existence of those bound states at the large m^i is also confirmed by the sign of 
scattering length aq in each channel [18, 20] based on the discussion in Refs. [15, 16]. 

In the spin-triplet ^Si channel, the bound state nature found in lattice calculations is consistent 
with the presence of deuteron in nature. The magnitude of binding energy, however, is significantly 
larger, and does not show a trend of decrease toward the experimental value as is taken lighter. 
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at least to = 300 MeV [21]. The trend is similar in the spin-singlet ^Sq channel. There is a 
bound state down to = 300 MeV, and it is not clear if or how this bound state disappears as one 
approaches the physical pion mass. 

A recent result for 8{p) in the 'Sq channel obtained by the HALQCD method [68] is shown 
in the right panel of Fig 5. If there is a bound state, 5{p) should start at 180 degrees at p = 0. The 
figure therefore shows that there is no bound state in this channel at = 410-700 MeV. In contrast 
to the results from Luscher’s method described above, this is consistent with experiment. Quantita¬ 
tively, however, the magnitude of phase shift in experiment is much larger than the values in Fig. 5. 
It would be important to perform a calculation at the physical to reproduce the experimental re¬ 
sult. Bound state is also not observed in the ^Si channel by the HALQCD method [73]. The 8{p), 
though increasing with decreasing nijt, does not support a bound state down to 470 MeV in 
Nf = 3 QCD. This situation is opposite to that with the Luscher’s method. 

We have to conclude that at present the results from the two methods in the NN channels do 
not agree, even for the presence or absence of bound states. 


AE(^Sj)[GeV] AE('S(,)[GeV] 



m^[GeV^] m^[GeV^] 


Figure 9: Energy shift in NN ^Si (left) and ^So (right) channels. Experimental value of deuteron is also 
plotted in left panel. Voo and Vmax express results in the infinite volume and on the largest volume, respec¬ 
tively. 


5.5 Systematic uncertainty 

We have seen that the results from the two methods agree with each other at least qualitatively 
for the 1 = 2 two-pion and //-dibaryon channels, while they differ in an essential way for the NN 
channels. A possible reason for this difference may be that the first two cases are simple two-body 
systems governed by potentials which are monotone functions of r. V (r) in the 1 = 2 two-pion 
channel decreases as r increases (Fig. 6), and the one in the //-dibaryon channel increases as r 
(Fig. 8). The NN potential, however, has a more complex structure (Fig. 5) with a repulsive core 
near the origin and an attractive tail at large distances. Thus, various systematic uncertainties 
might affect a balance toward bound state formation in a serious way. Let us then discuss possible 
systematic uncertainties in the two methods for the NN channels. 
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5.5.1 Liischer’s method 

For bound states, Liischer’s method is nothing but the traditional calculation of evaluating the 
mass of the ground state of hadrons from two-point correlation function, and extrapolating it to the 
infinite volume. Thus, the crucial point we need to worry about is how reliably the ground state 
energy is extracted for a given volume. Figure 10 presents the effective energy shift for the NN 
state from in the ^Si channel by Yamazaki et al. [19] and NPLQCD Collaboration [20]. While 
a clear plateau extends only over 3 or 4 time slices, the negative value of the energy shift appears 
quite convincing. 

One should note that the results above are obtained with only one type of source operator. It 
is a legitimate question to ask if one obtains the same result, when a different source operator is 
employed. Similarly, the dependence on the choice of the sink operator is also a potential issue. 
These issues have not been thoroughly investigated. Application of the variational method [46] 
would be worthwhile to settle these issues. 

There are two more sources of possible errors which, however, are more dynamical in nature. 
One is the dependence on the quark mass. It may well be that the bound state found in the ^ So 
channel for the pion mass down to nij^ k, 300 MeV [21] disappears at a smaller pion mass but 
heavier than the physical mass. 

Another is the dependence on the lattice spacing. It is possible in principle that the pion mass 
below which a bound state in the spin-singlet 'Sq channel disappears varies with the lattice 
spacing. It may require a sufficiently small lattice spacing before becomes non-zero or larger 
than the physical value. 

These possibilities can only be settled by repeating the calculations toward the physical pion 
mass and for smaller lattice spacings. 



Figure 10: Effective energy shift calculated by Yamazaki et al. [19] (left) and NPLQCD Collaboration [20] 
(right) in NN ^Si channel. 


5.5.2 HALQCD method 

The HALQCD method relies crucially on the potential V (r) extracted from the NBS wave 
function. Hence possible systematic errors center around the meaning of the potential and how the 
final physics resulls may or may nol depend on fhe ambiguifies in fhe pofenfial. Lef us discuss 
possible sources of syslemalic errors and ambiguifies one by one. 
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(i) Sink operator dependence of V(r) and d(p) 

The potential V (r) will in general depend on the choice of the sink operator of CNN{r,t) in eq.(5.1). 
In SU(2) gauge theory, Takahashi et al. [79] examined this issue, and found that V (r) in the small 
r region largely depends on the choice of the sink operator, while the dependence is small for large 
r. To what extent such a dependence affects the phase shift extracted from V (r) was not tested in 
this study, but this is certainly a point which requires a detailed check. 

(ii) Velocity expansion of U{f,r') 

If higher order of velocity expansion is not negligible in eq.(5.4), V{f) obtained from eq.(5.3) 
or (5.5) depends on the relative momentum p of the two-nucleon state. The dependency was 
investigated by Murano et al. [80] at p ~ 0 and 250 MeV, and it was concluded that the V(r) 
obtained from two values of p agree within the statistical error. 

(iii) Volume dependence of V (r) 

If the extracted potential V(r) has finite volume effects, it may affect S(p) and hence the presence 
or absence of bound states. Thus, finite volume dependence was investigated by HALQCD Col¬ 
laboration [72] at = 1 GeV in the spatial extent of 2-3 fm, and it was concluded that the finite 
volume effects in V(r) is negligible within the statistical error. 

(iv) Lattice spacing dependence of V (r) 

Lattice spacing dependence of V(r) was examined at = 1.1 GeV for three different lattice 
spacings, a = 0.11-0.22 fm, by HALQCD Collaboration [81]. They found that lattice spacing 
dependence is small in the larger r region, but it is large in the small r region, r;^0.5 fm, and then 
S(p)’s obtained at different lattice spacings differ at large p region. 

(v) Determination of mf/ 

If the value of contains systematic errors, which is even a few MeV, it corresponds to a constant 
shift of V (r) for all values of r, since 


V(r) 


(S- 

CNN{r,t) 




(5.6) 


The tiny constant shift of V{r) might be crucial in some case, e.g., identification of existence of 
bound state. 

(vi) V(r) in large r region 

It is hard to extract V (r) in the large r region precisely, because the noise-to-signal ratio of V (r) 
increases with r. Since the large r region is important to determine how long attractive tail of the 
potential continues, this uncertainty might affect existence of bound state as pointed out by the 
plenary speaker in the last conference [82]. 

To sum up the considerations above, most uncertain with the HALQCD method seems to be 
the ambiguities in the potential for small values of r as written in (i), (iii), and (iv). The small r 
region is not relevant to calculate 5{p) at small p, while it becomes important to calculate 5{p) 
at large p, and also to test the existence of a bound state. Thus, those uncertainties in the small r 
region might cause the difference of the result obtained from the two methods. 


6. Light nuclei 

A first attempt to simulate the formation of helium nuclei in lattice QCD was reported in 
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Ref. [83]. In this work, the binding energies of the ^He and ^He nuclei were calculated in quenched 
QCD at = 800 MeV, by examining the volume dependence of the energy shift. A serious 
computational problem with nuclei calculations is huge number of Wick contraction. For helium 
nuclei [83], this problem was overcome by omitting calculations of redundant contractions under 
symmetries of interpolating operator, and by utilizing blocks of three quark propagators. Recently, 
more efficient calculation methods [84, 85, 86] have been proposed. In this conference a prelimi¬ 
nary result in this direction [87] is reported. 

The exploratory study of helium nuclei has been followed by calculations in A/ = 3 QCD at 
niT, = 810 MeV [20], and 2 + 1 QCD at = 510 MeV [19] and = 300 MeV [21]. Results 
for the binding energies for ^He and ^He nuclei are summarized in Fig. 11. For both the nuclei, 
the quenched and Nf = 3 calculations at ~ 800 MeV give different results. While it might be a 
dynamical quark effect, precise understanding is lacking at present. At m„ = 300, the lightest pion 
explored so far, the binding energy of ‘^He is roughly consistent with the experiment, while that 
for ^He is about three times larger than the experiment. The discrepancy between the lattice and 
experimental results might be caused by a large in the calculation. Future calculations closer to 
the physical nijc should tell if this expectation is justified. 

In fhe righf panel of Fig. 11 fhe resulf from HALQCD mefhod for '^He in Nf = 3 QCD [73], 
giving a very small binding energy, is also plotted. Resulfs for ^He is nof available from fhe 
HALQCD mefhod. 
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Figure 11: Binding energies for ^He (left) and ^He (right) channels. Experimental values are also plotted. 
Voo and Vmax express results in the infinite volume and on the largest volume, respectively. 


Let us list some related studies for nuclei. NPLQCD Collaboration calculates magnetic mo¬ 
ments for light nuclei [88] and binding energies for quarkonium-nucleus bound states [89]. HALQCD 
Collaboration calculates the binding energy of spin-2 NQ. bound state [90] and medium-heavy 
nuclei [73] from potentials obtained by the HALQCD method. The nuclei calculation has been 
extended to theories beyond the standard model. Detmold et al. [91] calculate nuclei of two to 
four particles in two-flavor SU(2) gauge fheory fo explore fhe imporfance of nuclear physics in fhe 
sfrongly-inleracfing dark matter models. 
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7. Conclusion 

In this report we review reeent works for seattering and deeay ehannels. The seattering length 
in 7 = 2 two-pion ehannel is preeisely determined using ehiral extrapolation with ChPT formula. 
Still eareful estimation of systematie errors are neeessary, some reeent results have smaller error 
than the experimental results. Thus, the ealeulation in this ehannel enters preeision measurement 
era as in hadron masses and deeay eonstants. The ealeulation of seattering length in 7 = 1/2 K% 
ehannel is more diffieult than the one in 7 = 2 two-pion ehannel, so that only a few ealeulations 
are earried out. The resonanee mass and deeay width are extraeted from seattering phase shift. In 
7=1 two-pion ehannel, the p meson mass and its eoupling to two pions are ealeulated. Both the 
physieal quantities are roughly eonsistent with the experiment within large error of lattiee results. 

Comparisons of the results from Liiseher’s method and HALQCD method are diseussed. In 
7 = 2 two-pion ehannel, the results from the two methods agree with eaeh other numerieally. While 
the binding energies of 77-dibaryon differ in larger m^i region, the two method give the eonsistent 
eonelusion that the 77-dibaryon exists in larger niji than physieal one. However, in the NN ehannels, 
the two methods obtain opposite results in the point of view of the presenee of bound states. Pos¬ 
sible uneertainties in eaeh method are diseussed in the NN ehannels. Contaminations from higher 
exeited states, larger quark mass, and finite lattiee spaeing eould affeet the result from Liiseher’s 
method. For HALQCD method, several possible uneertainties are listed. Uneertainties of the po¬ 
tential in small r region might explain why the bound state is not formed in this method. In order 
to understand the differenee between the results obtained from Liiseher’s method and HALQCD 
method, it would be neeessary to earry out more detail investigations for uneertainties in both the 
methods. Furthermore, understanding of relation between the potential evaluated from NBS wave 
funetion, in other words, Fourier transformation of off-shell seattering amplitude, and phase shift 
in quantum field fheory mighf shed lighf on fhe differenee. 

The helium nuelei are ealeulafed in several works using Luseher’s mefhod. In fhe ^He ehannel 
fhe binding energy is roughly eonsisfenf wifh fhe experimenf, while in fhe ^He ehannel if is larger 
fhan fhe experimenf. The diserepaney of fhe binding energy from fhe experimenf is eonsidered 
fo be mainly eaused by large in fhe ealeulation. Therefore, ealeulafions af smaller ruji < 300 
MeV would be neeessary fo sfudy whefher lattiee ealeulafion reproduees fhe experimenfal binding 
energy or nof. 
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